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ABSTRACT 



Aims. We present deep H.E.S.S. observations of the supernova remnant (SNR) RX J1713. 7-3946. Combining data of tiiree years - from 2003 
to 2005 - we obtain significantly increased statistics and energy coverage as compared to earlier 2003 & 2004 results. 
Methods. The data are analysed separately for the different years. 

Results. Very good agreement of the gamma-ray morphology and the differential spectra is found when comparing the three years. The com- 
bined gamma-ray image of the 2004 and 2005 data reveals the morphology of RX J1713. 7-3946 with unprecedented precision. An angular 
resolution of 0.06° is achieved, revealing the detailed structure of the remnant. The combined spectrum of all three years extends over three 
orders of magnitude, with significant gamma-ray emission approaching 100 TeV. The cumulative significance above 30 TeV is 4.8o", while for 
energies between 113 and 294 TeV an upper limit on the gamma-ray flux of 1.6 x 10"'* cm"^ s"' is obtained. 

Conclusions. The energy coverage of the H.E.S.S. data is presumably at the limit of present generation Cherenkov telescopes. The measure- 
ment of significant gamma-ray emission beyond 30 TeV formally implies the existence of primary particles of at least that energy. However, for 
realistic scenarios of very-high-energy gamma-ray production, the Inverse Compton scattering of very-high-energy electrons and n° decay fol- 
lowing inelastic proton-proton interactions, the measured gamma-ray energies imply that efficient acceleration of primary particles to energies 
exceeding 100 TeV is taking place in the shell of the SNR RX J1713.7-3946. 

Key words, acceleration of particles - ISM: cosmic rays - gamma rays: observations - ISM: supernova remnants - gamma rays: individual 
objects: RX J1713.7-3946 (G347.3-0.5) 
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These observations had wobble offsets of 0.5° in Right Ascension and Declination, less than the radius of the SNR, and are 



summarised in one row. 

Table 1. Summary of H.E.S.S. observations of RX J1713. 7-3946 conducted during three years. For each year, the targeted po- 
sition is given together with the wobble offsets in Right Ascension (ffjaooo) and DecHnation ((5j2ooo)- Adding the wobble offsets 
to the target coordinates, the actual pointing position is obtained. For each pointing position, the dead-time corrected observation 
time {Live time) is given in hours. Data recorded under bad weather conditions are excluded. The columns Data set I-III sum- 
marise observation times of data sub-sets used throughout the paper. Data set I is used for Fig. |2l to obtain optimum angular 
resolution, the 2003 data are disregarded and a zenith-angle cut at 60° is applied (the latter is only relevant for the 2005 data). 
Data set II is used for spectral comparisons of the different years (cf. FigO). Observations with wobble offsets of 0° in 2004 and 
2005 are disregarded for this purpose. Data set III comprises all data, and is used for the generation of the combined spectrum 
(cf. Fig.©. 



1. Introduction 

The energy spectrum of cosmic rays measured at Earth exhibits 
a power-law dependence over a broad energy range. Starting at 
a few GeV (1 GeV =10^ eV) it continues to energies of at least 
10^" eV. The power-law index of the spectrum changes at two 
characteristics energies: in the region around 3 x 10'^ eV - the 
knee region - the spectrum steepens, and at energies beyond 
10^^ eV it hardens again. This latter feature is known as the an- 
kle. Up to the knee, cosmic rays are believed to be of Galactic 
origin, accelerated in shell-type supernova remnants (SNRs) - 
expanding shock waves initiated by supernova explosions (for 
a recent review see Hillas (2005.) ). However, the experimental 
confirmation of an SNR origin of Galactic cosmic rays is dif- 
ficult due to the propagation effects of charged particles in the 
interstellar medium. The most promising way of proving the 
existence of high-energy particles in SNR shells is the detec- 
tion of very-high-energy (VHE) gamma rays {E > 100 GeV), 
produce d in interactions of cosmic rays close to their accelera- 
tion site dOrurv et alJll994 . 

Recently H.E.S.S. - a VHE gamma-ray instrument 
consisting of four Imaging Atmospheric Cherenkov 
Telescopes - has detected two shell -t ype S NRs, 
RX J1713.7-3946 ("Aharonian et al.' '2004b', 2006b) and 
RX 10852.0^622 (Aharonian et al . 2005a). The two objects 
show an extended morphology and exhibit a shell structure, 
as expected from the notion of particle acceleration in the ex- 
panding shock fronts. Both objects reveal gamma-ray spectra 
that can be described by a hard power law (with photon index 
r ~ 2.0) over a broad energy range. For RX J1713. 7-3946 
significant deviations from a pure power law at larger energies 



are measured ( Aharonian et al. 2006b'). While it is difficult 
to attribute the measured VHE gamma rays unambiguously 
to nucleonic cosmic rays (rather than to cosmic electrons, 
which would certainly also be accelerated in the shock front), 
the measured spectral shapes favour indeed in both cases 
- for RX J1713.7-3946 and RX J0852.0-4622 - a nucle- 



onic cosmic -ray origin of the gamma rays (lAharonian et al 



|2006a,b). In the case of RX J171 3.7-3946 in additio n a narrow 
shock filament seen in X-rays (iHiraga et alj l2005h indicates 
strong amplification of the magnetic field at least in one region 
of the rim ( Bere zhko & Volk 2006 ). If such an amplified 
magnetic field exists throughout the main volume of the 
SNR - the region for which VHE gamma-ray data is presented 
here - and if consequently high magnetic field values are found 
not only in one shock filament, but on a large part of the shock 
surface, a leptonic origin of the VHE gamma rays becomes 
increasingly unlikely just based on the absolu te level of X-ray 



and ga mma-ray flux of RX J1713. 7-3946 dAharonian et al 



2006b) 



Apart from the first unambiguous proof of multi-TeV parti- 
cle acceleration in SNRs ( Aharonian et al, 2006a,b), the ques- 
tion of the highest observed energies remains an important one. 
Only the detection of gamma rays with energies of 100 TeV 
and beyond provides experimental proof of acceleration of 
primary particles, protons or electrons, to even higher ener- 
gies of 1 PeV and beyond. T he spectrum of th e whole SNR 



RX J1713. 7-3946 reported in lAharonian et al 



prises data of the 2004 observation campaign of H.E.S.S.. It 
extends to energies of 40 TeV. Here we present a combined 
analysis of H.E.S.S. data of RX J1713.7-3946 recorded in 
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three years, in 2003 during the construction and commission- 
ing phase of the system, and in 2004 and 2005 with the full 
H.E.S.S. array. A comparison of the three data sets demon- 
strates the expected steady emission of the source as well as 
the stability of the system during the first three years of run- 
ning. Special emphasis is then devoted to the high-energy end 
of the combined spectrum. 

2. H.E.S.S. observations 

The High Energy Stereoscopic System (H.E.S.S.) consists of 
four identical Cherenkov telesco pes that are ope rated in the 



tour Identical CherenRov telesco pes that are ope rated m the 
Khoma s Highland of Namibia (iHofmannI 2005 ). The tele 
scopes ( Bernlohr et al. 2003 : Cornils et al. j2003h are 13 m in 
diameter, each with a mirror area of 107 m^. During nor- 
mal operation the system is run in a coincidence mode which 
requires a trigger from at least two out of the four tele- 
scopes ( Funk et al. 112004*). The cameras consist of 960 photo 



multiplier pixels and cover a 5° field of view (IVincent et al 
2003h . The resulting FWHM « 4° of the system field-of-view 



response makes H.E.S.S. the currently best suited experiment 
in the field for the study of extended VHE gamma-ray sources 
such as young Galactic SNRs. At zenith, the energy thresh- 
old is about 100 GeV and for point sources an energy resolu- 
tion of 15% is achieved. The angular resolution for individual 
gamma rays is better than 0.1° and the point source sensitivity 
reaches 1 % of the flux of the Crab nebula for long exposures 
(~ 25 hours). 

The H.E.S.S. observation campaign of RX J1713.7-3946 
started in 2003. The data were recorded between May and 
August 2003 during two phases of the commissioning of the 
telescope system. During the first phase, two telescopes were 
operated independently with stereoscopic event selection done 
offline using GPS time stamps to identify coincident events. 
During the second phase, also using two telescopes, coinci- 
dent ev ents were selected in hardware using the array level 
trigger dFunk et al.ii2004) . The observations were performed in 
Declination wobble mode around the northwest shell of the 
SNR, the alternating wobble offset in Declination was 0.5°. 
The zenith angle of observations varied from 15° to 30° with 
a mean of 24°. The analysis of this first data set revealed ex- 
tended gamma-ray emission resembling a shell structure, very 
similar to the X-ray image. It was actually the first ever re- 
solved imag e of an astronomical so urce obtained with VHE 
gamma rays dAharonian et al. 2004bV The spectrum was well 
described by a hard power law with energies from 1 to 10 TeV. 

In 2004, observations were conducted with the full tele- 
scope array. From April to May, most of the data were recorded 
in wobble mode, this time around the SNR centre with an off- 
set of 0.7° in Right Ascension and Declination aiming at more 
uniform coverage of the whole SNR and, important for anal- 
ysis purposes, fully encompassing the SNR with the four ob- 
servation positions. The zenith angle of observations ranged 
from 16° to 56° with a mean of 26°. The H.E.S.S. data en- 
abled analysis of the gamma-ray morphology and the spectrum 
of the remnant with unprecedented precision ( Aharonian et al.l 
12006b'). A very good correlation was found between the X-ray 
and the gamma-ray image. The differential spectrum was mea- 



sured from 200 GeV up to 40 TeV. A deviation from a pure 
power law was found at high energies. A spatially resolved 
spectral study revealed no significant changes of spectral shape 
across the SNR despite flux variations by more than a factor of 
two. 

The 2005 observation campaign was aiming at extending 
the energy coverage of the spectrum to as high energies as pos- 
sible. Therefore the observations, carried out from beginning 
of September to November, were preferentially pursued at large 
zenith angles, up to values of 70°, to make use of the drastically 
increased effective collection area of the experiment at high en- 
ergies. The mean zenith angle of observations was 51°. As in 
2004, RX J1713.7-3946 was observed in wobble mode with 
an offset of 0.7° in Declination and Right Ascension. Analysis 
of these data are for the first time presented in the following. 
A summary of the observations conducted during three years 
with H.E.S.S. is given in Table[T] 

3. Data Analysis 

The RX J1713. 7-3946 data presented here are cali- 
brat ed according to the sta ndard H.E.S.S. calibration meth- 
ods ( Aharonian et al. 2004al) . For the backgro und suppression, 
cuts on scaled image parameters are applied ( Aharonian et alJ 
i2005b ). The shower reconstruction is based on image param- 
eters (Hillas parameters) and corresponds, un less otherwise 
stated, to algorithm 1 of iHofmann et al. (Il999h : the intersec- 
tion point of the image axes in a common camera coordinate 
system yields the shower impact position on ground and the di- 
rection of the primary. A cut on the minimum size of camera 
images is applied to assure that only well defined images are 
included in the analysis. For the 2003 two-telescope data, the 
cut is applied at a rather large value of 300 photo-electrons. In 
the commissioning phase of the experiment, this served to dra- 
matically reduce the number of background events, but it also 
homogenises the whole data set, which was recorded with two 
different hardware configurations, thereby reducing systematic 
uncertainties. Moreover, the angular resolution improves when 
including only well defined images in the an alysis. The 2004 



and 200 5 data are analysed as discussed in lAharonian et al 



(l2006bl) . "or spectral analysis, a loose cut on the minimum im- 
age size at 80 photo-electrons is applied. For studies of the 
gamma-ray morphology, the cut is increased to 200 photo- 
electrons yielding superior angular resolution of the order of 
0.08° and better background suppression. 

For the subtraction of the irreducible cosmic-ray back- 
ground, separate approaches are taken for the generation of 
gamma-ray spectra and images. The preferred background- 
estimati on method for spec tral analysis is the reflected- region 
model (IHinton et alJ l2005h . The background estimate is de- 
rived from a region of the same size and shape as the source 
region, reflected at the system pointing direction. To assure 
non-overlapping source and background-control regions, this 
approach can only be applied if the observation positions have 
been chosen outside the nominal gamma-ray source region. As 
can be seen from Table [T] this is not true for the whole 2003 
data and parts of the 2004 and 2005 observations. For these 
data, an ON/OFF-hackgmund model is applied instead. From 
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Year 




^68 


ON 


OFF 


a 


Significance (cr) 


Live time (hrs) 


2003 


24° 


0.083° 


3194 


1764 


1.00 


21 


18.1 


2004 


27° 


0.075° 


18728 


11039 


1.05 


41 


33.1 


2005 


44° 


0.082° 


10277 


5124 


1.15 


33 


29.7 



Table 2. Summarised are the event statistics of the whole SNR and corresponding angular resolutions for the years 2003, 2004, 
and 2005. The 2004 and 2005 data sets correspond to Data set I of Table[T] The average zenith angle (0z) is determined from all 
events reconstructed in the SNR region. For the angular resolution {R(,s), the 68% containment radius of the simulated gamma-ray 
point-spread function, matched to the particular data set, is used as figure of merit. The other columns give the number of signal 
events from the SNR region (ON), the number of background events (OFF), the normalisation factor between ON and OFF counts 
(a), and the corresponding significance and live observation time, a is in general defined as the ratio of the effective exposure 
integrated in time and angular space of the ON and OFF region. Note that the analysis of the 2003 data is adopted to match the 
system configuration of this year The nominal analysis is applied for the 2004 and 2005 data. For 2005, only data recorded at 
zenith angles less than 60° are included (therefore the mean zenith angle decreases). The event statistics are determined with the 
ON/OFF approach for 2003 and the reflected-region method for 2004 and 2005. In the latter two years, also ON runs with wobble 
offsets < 0.7° (cf. Table[TJ are included and hence a > 1. 



the complete set of H.E.S.S. observations without gamma-ray 
signal, OFF runs for background estimation are selected with 
zenith-angle distributions matching that of the ON runs as close 
as possible. 

For im age generation, the fi e/(i-of-v/ew-background model 
is applied (iHinton et alJ 120051) . It models the background by 
means of a system acceptance model determined from the full 
set of H.E.S.S. OFF runs. The normalisation is calculated using 
the whole field of view excluding regions of known gamma- 
ray emission. Note that the background-subtracted gamma- 
ray images shown throughout this paper are smoothed with a 
Gaussian to reduce statistical fluctuations. The resulting im- 
ages are in units of gamma-ray excess counts per Gaussian 
sigma. They are corrected for the falloff of the system accep- 
tance towards the edges of the field of view which results from 
a smaller detection efficiency far from the pointing centre. 

When determining spectra of the whole SNR, a circu- 
lar region of 0.65° radius is used here, centred at amm - 
17hl3m33.6s, 5mm = -39°45'36". 



4. Gamma-ray Morphology 

The gamma-ray morphology as obtained from the H.E.S.S. 
data of three years is seen in the upper panel of Fig. [T] For the 
2003 data, a special set of analysis cuts was applied adopted 
to the two-telescope data (cf. Sec.[3]l. For the 2005 image, only 
observations at zenith angles smaller than 60° are taken into ac- 
count (reducing the available observation time by ~ 10 hours, 
cf. TablefUl. For larger zenith-angle observations, the geometri- 
cal reconstruction worsens, deteriorating the angular resolution 
of the resulting image. 

The images of Fig.[T]are readily comparable. They are cor- 
rected for system acceptance, which is different for the dif- 
ferent data because of the zenith-angle dependence of the ac- 
ceptance and the intermediate system configuration of 2003. 
Very similar angular resolutions are achieved for all years, see 
Table |2l where relevant parameters are listed. From the visual 
impression the three images shown in the figure are very sim- 
ilar Within statistics, good agreement is achieved, as can be 
seen from the one-dimensional distributions shown in the lower 



panel of Fig.[Tl which have the advantage that statistical errors 
on the measurement can be taken into account for the com- 
parison. Shown from left to right are a slice along a thick box 
(shown in the upper panel of Fig. [TJ, an azimuthal profile of 
the shell region, and a radial profile. All the distributions are 
generated from the non-smoothed, acceptance-corrected excess 
images, very finely binned such that binning effects are negli- 
gible. Clearly, there is no sign of disagreement or variability, 
the H.E.S.S. data of three years are well compatible with each 
other. 

The combined H.E.S.S. image is shown in Fig. [2] Data 
of 2004 and 2005 are used for this smoothed, acceptance- 
corrected gamma-ray excess image {Data set I in Table [T]|. In 
order to obtain optimum angular resolution, a special analysis 
is applied here. In addition to the image-size cut of 200 photo- 
electrons, the cut on the minimum event multiplicity is raised 
to three telescopes (disregarding the 2003 data for this purpose 
completely). Moreover, an advanced reconstruction method is 
chosen. It takes Monte-Carlo error estimates on image parame 



ters into acc ount an d is based on algorithm 3 of iHofmann et al 



fl999) (see iBergd i2006) for studies of this analysis tech 



nique). The image corresponds to 62.7 hours of dead-time 
corrected observation time. With the reflected-region method, 
12961 ON events from the region associated with the SNR are 
accumulated, and 5710 OFF events (normalisation a - 1. 10). 
Hence, 6702 gamma-ray excess events are measured with a 
statistical significance of 48cr. An angular resolution of 0.06° 
(3.6') is achieved. For comparison, the resolution obtained with 
the standard geometrical reconstruction method and a three- 
telescope multiplicity is 0.07° with similar event statistics. 
With a two-telescope multiplicity cut, the resolution with the 
standard reconstruction is 0.08° (with 28879 ON, 16070 OFF 
events, cc = 1.1, and a significance of 53 cr). 

The image in Fig.|2lconfirms nicely the pub Ushed H.E.S.S. 
measurements jAharonian et al ] |2004bl.l2006 b'), with 20% bet- 
ter angular resolution and increased statistics. The shell of 
RX J1713. 7-3946, somewhat thick and asymmetric, is clearly 
visible and almost closed. As can be seen from the left-hand 



' Note that pure ON runs with wobble offsets < 0.7° are included 
in Data set I (cf . Table[TJ and hence a > I. 
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Fig. 1. Upper panel: H.E.S.S. gamma-ray excess images from the region around RX J1713. 7-3946 are shown for each year 
separately for comparison. From left to right, images are generated from data of 2003, 2004 and 2005. The images are corrected 
for the decline of the system acceptance with increasing distance to the SNR centre. All three images are smoothed with a 
Gaussian of 2', the Unear colour scale is in units of excess counts per smoothing radius. The dashed box (dimensions 2° x 0.6°) 
and ring (ri = 0.3°, r2 - 0.5°) are used for obtaining the one-dimensional distributions shown in the lower panel. Lower panel: 
One-dimensional distributions for the data of three years, all generated from the raw, non-smoothed and acceptance-corrected 
gamma-ray excess images. Histograms of 2004 and 2005 are scaled to the area of the 2003 histogram to account for differences 
in the event statistics. Left: Slices taken within a rotated box running through the SNR region. Plotted are events versus angular 
distance to the centre, projected onto the axis running through the SNR centre, rotated by 45° anti-clockwise with respect to the 
RA axis. Middle: Azimuth profiles integrated in a thick ring covering the shell of RX J1713. 7-3946. The azimuthal angle of the 
events is calculated with respect to the SNR centre (Qrj2ooo = 17hl3m33.6s, 5j2()oo - -39°45'36")- 0° corresponds to the west 
part of the shell, 90° is north or upward, -90° is south or downward. Right: Radial profiles around the centre of the SNR. Plotted 
are excess events per unit solid angle as a function of the distance r to the SNR centre. 



side of the figure, when integrating signal and background 
events in a circle of 0.1° radius around each trial point-source 
position, significant gamma-ray emission is found throughout 
the whole remnant. Even in the seemingly void south-eastern 
region it exceeds a level of 8 standard deviations. The gamma- 
ray brightest parts are located in the north and west of the SNR. 
The similarity of gamma-ray and X-ray morphology, which 
was already investigated in detail in Aharonian et al. (2006b) 
for the 2004 H.E.S.S. data, is again demonstrated on the right- 
hand side of Fig.|2] where ASCA X-ray contours are overlaid 
on the H.E.S.S. image. 

5. Gamma-ray Spectrum 

The gamma-ray spectra measured with H.E.S.S. in three con- 
secutive years are compared to each other in Fig. [3] The 2003 
spectrum is obtained from an ON/OFF analysis, with the set 



of special two-telescope cuts mentioned above. Note that these 
cu ts were also app lied to obtain the spectrum shown in Fig. 3 
of lAharonian et aL I (l2004bh . which stops at 10 TeV. Here, how- 
ever, the 2003 spectrum extends to energies beyond 30 TeV. 
The difference between the two analyses is the energy range of 
simulations used to generate effective gamma-ray detection ar- 
eas (needed for spectral analysis). In the old analysis, gamma 
rays were simulated up to 20 TeV, permitting energy recon- 
struction only up to ~ 10 TeV (allowing for a maximum re- 
construction bias of 10%). Here, in the present analysis, sim- 
ulations up to 100 TeV are available for zenith angles smaller 
than 60°, up to 200 TeV for angles from 60° to 63°, and up to 
400 TeV for zenith angles up to a maximum of 70°. Hence the 
increased energy coverage. Note that good agreement is found 
between the 2003 spectrum shown he re and the one published 
previously in lAharonian et al.l (l2004bl) in the energy range from 
1 to 10 TeV. 
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Fig. 2. Two versions of the combined H.E.S.S. image from the 2004 and 2005 data. Shown is in both cases an acceptance- 
corrected gamma-ray excess image. The images are smoothed with a Gaussian of 2'. A simulated point source as it would appear 
in this data set is shown in the lower left-hand corner of both images (labeled PSF). It is smoothed with the same Gaussian, the 
cr of 2' is denoted as black circle in the insets. The linear colour scale is in units of excess counts per smoothing radius. Note that 
for the 2005 data, only data recorded at zenith angles less than 60° are taken into account. On the left-hand side, the overlaid 
light-gray contours illustrate the significance of the different features. The levels are at 8, 18, and 24cr. The significance at each 
position has been calculated for a point-source scenario, integrating event s in a circle of 0. 1 ° rad ius around that position. On the 
right-Iiand side, ASCA contours are drawn as black lines (1-3 keV, from lUchivama et al.l (120021) ) for comparison. 



Year 


ON 


OFF 


a 


Significance (cr) 


Live time (hrs) 


2003 


3194 


1764 


1.0 


20.5 


18.1 


2004 


107494 


93906 


1.0 


30.3 


30.6 


2005 


71276 


60175 


1.0 


30.6 


36.8 



Table 3. Comparison of event statistics from the SNR region from three years of data. The numbers result from the spectral 
analysis of Data set II (cf. Table [TJ, shown in Fig. [3] Given are the number of signal (ON) and background (OFF) counts, the 
normalisation factor a, the statistical significance of the gamma-ray excess (cr) and the observation time. For the 2003 data, 
the special two-telescope analysis with a cut on the minimum size of camera images at 300 photo-electrons was applied. The 
background estimate in this case is derived with the ON/OFF analysis. For 2004 and 2005, the nominal spectral analysis with a 
cut at 80 photo-electrons was used together with reflected-region background model. 



The spectra determined from the 2004 and the 2005 data 
in Fig. [3] are obtained with the reflected- region-hackground 
model. Therefore, data where the observation position was 
within the SNR region are disregarded. For the purpose of 
comparison of the different data sets this approach seems rea- 
sonable, no attempt to analyse the remaining data with an 
ON/OFF -background approach is pursued. The corresponding 
event statistics for the spectra shown in the figure are listed in 
Tabled 

In order to compare data recorded in different years, a cor- 
rection for the variation of optical efficiency of the telescope 
system must be applied. The efficiency degrades with time, 
mostly due to degradations of miiTo r reflectivity. As described 
in detail in lAharonian et al.l (l2006cl) . this worsening of the ac- 
tual efficiency with respect to the simulated one causes a shift 
in the absolute energy scale. This shift can be corrected us- 



ing measured images of local muons, for which the light yield 
is predictable. Based on the prediction and the simulated light 
yield, an average energy correction factor is determined for the 
data of each of the three years separately. The resulting av- 
erage values are 1.12 for 2003 and 2004, and 1.30 for 2005. 
These correction factors are used to correct the reconstructed 
energies thereby enabling direct comparisons between differ- 
ent years. Note that a correction factor is needed already for 
the first data set of 2003 since the Monte-Carlo simulations re- 
fer to new mirrors, but in 2003 the first H.E.S.S. telescope was 
already one year old. In 2004, the total optical efficiency of the 
system remained the same because of the inclusion of two tele- 
scopes with nominal efficiency, thereby cancelling the aging 
effects of the first two telescopes. 

The spectra shown in Fig.|3]are after correction. Very good 
agreement is found between the different years. The measured 
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Fig. 3. Comparison of H.E.S.S. spectra from the years 2003, 
2004, and 2005 {Data set II, Table [B- The three spectra are 
shown in an energy-flux representation - flux points have been 
multipUed by E^. The black curve is shown for reference. It 
is the best fit of a power law with exponential cutoff" to the 
combined data, where the cutoff is taken to the power of [3 - 
0.5: dNIdE = /(, E'^ exp [-{EjE^) P=^-^). Note that flux points 
are corrected for the degradation of the optical efficiency of the 
system. The energy threshold of ~ 1 TeV in the 2003 data is 
due to the two-telescope operation mode and the application of 
a stringent cut on the minimum camera image size. 



spectral shape remains unchanged over time. The absolute flux 
levels are well within the systematic uncertainty of 20%. As ex- 
pected for an object like RX J1713. 7-3946, no flux variation 
is seen on yearly timescales. Clearly, the performance of the 
telescope system is under good control, the correction of the 
optical degradation by means of energy correction factors de- 
ter mined from "muon efficie ncies" works reasonably well (see 
also lAharonian et al.l (l2006cl) ). Note that without correction of 
aging effects, flux differences between 2004 and 2005 are on 
the order 40%. 

The combined data of three years are shown in Fig. |4] 
This energy spectrum of the whole SNR region corresponds 
to 91 hours of H.E.S.S. observations {Data set III, Table [T]i. 
It is generated analysing the 2003 data separately, with the 
ON/OFF approach and the two-telescope analysis. The 2004 
and 2005 data are analysed together, with the reflected-region 
background and the nominal 80-photo-electrons cut. As shown 
in Table [TJ a fraction of the data was recorded with wobble 
off'sets smaller than 0.7°. For this part, the ON/OFF method is 
applied. Average energy-correction factors for each of the three 
subsets of data separately are determined as explained above. 
Having analysed the data separately to obtain suitable back- 
ground estimates for the SNR region, spectra are then com- 
bined to yield the final spectrum shown in Fig. |4] Systematic 
checks included the application of tighter cuts on the image 
amplitude to the 2004 and 2005 data and separate analysis of 
data recorded under small and large zenith angles (below and 
above 50°). While the spectra determined with diff'erent cuts 



Fig. 4. Combined H.E.S.S. gamma-ray spectrum of 
RX J1713.7-3946 generated from data of 2003, 2004, 
and 2005 {Data set III, Table [TJ. Data are corrected for the 
variation of optical efficiency. Error bars are ±lcr statistical 
errors. These data might be described by a power law with ex- 
ponential cutoff of the form dN/dE = h)E^^ exp^-(£'/£'c)^)- 
The best fit result (black solid line) is given here for y6 = 0.5 
(fixed), r = 1.8, and E^ = 3.7 TeV (cf. Table |4]for the exact 
values). Note that the fit function extends as dashed black line 
beyond the fit range for illustration. For comparison, the best 
fit of a power law with exponenti al cutoff" and = 1 , obtained 
solely from the 2004 data (Aharonian et al.ll2006bl) . is shown 
as dashed red line. A model-independent upper limit, indicated 
by the black arrow, is determined in the energy range from 113 
to 300 TeV. 



are fully compatible, a slight ffux overestimation is found for 
the large zenith-angle data, on the 10% level. The investigations 
of systematic uncertainties at the largest zenith angles are still 
underway, but since the eff'ect on the final spectrum is small, 
< 5%, the combined spectrum given here includes all data, up 
to zenith angles of 70°. 

The energy binning of the differential ffux shown in the 
figure is chosen to be 12 bins per decade. For the final two 
points of the spectrum, beyond energies of 30 TeV, the binning 
is three times coarser, 4 bins per decade, accounting for de- 
creasing event statistics at the highest energies. For the actual 
positioning of the ffux points w'ithin an energy bin, the method 
proposed in Lafferty & WyattI (1995) is adopted: the point ap- 
pears at the energy value, where the ffux value predicted by 
an eff'ective-area weighted model spectral shape (a power law 
with exponential cutoff) is equal to the mean ffux value over the 
energy bin. Note that this is only relevant for wide bins. Here, 
for the spectrum of Fig. |4] the procedure results in ffux points 
that are placed within 1 % of the central energy value of the bin. 
Only for the two last points the difference is considerable, they 
end up at 7 and 12% smaller energy values than the bin centre. 

The combined H.E.S.S. spectrum of RX J1713.7-3946 
shown in Fig. |4] extends over almost three decades in en- 
ergy, and is compatible with previous H.E.S.S. measurements. 
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Table 4. Fit results for different spectral models. The fit range is chosen from 0.3 to 1 13 TeV. The differential flux normalisation 
Iq is given in units of 10"'^ cm"^ s"' TeV"'. Shown are a power-law model (row 1), a power law with an exponential cutoff 
(row 2, 3, 4; the cutoff energy is given in TeV), a power law with an energy dependent photon index (row 5), and a broken 
power law (row 6; in the formula, the parameter S - 0.6 describes the sharpness of the transition from Fi to F2 and it is fixed in 
the fit). Note that when fitting a broken power law to the data, some of the fit parameters are highly correlated. 



Fit Fomula 


Fit Parameters (d.o.f.) 




lo = 20.5 ±0.4 F = 2.32 ± 0.01 


145.6 (25) 


loE-^ exp (-(£/£,/) 


/o = 21.3 ±0.5 F = 2.04 ±0.04 £, = 17.9 ± 3.3 /3 = 1.0 
/o = 34.1 ±2.5 F= 1.79 ±0.06 = 3.7 ± 1.0 /? = 0.5 
/o= 40.5 ±1.5 F= 1.74 ±0.02 = 2.3 ± 0.2 /? = 0.45 


39.5 (24) 
34.3 (24) 
34.2 (24) 


/o S-r+Z' tag'; 


lo = 20.6 ±0.5 F = 2.02 ± 0.04 p = -0.29 ± 0.03 


38.8 (24) 


lo (E/Eb)-'' (i+(£/£b)"')'*""'" 


/o = 0.5 ±0.4 Fi= 2.00 ±0.05 F2 = 3.1±0.2 = 6.6 ± 2.2 


29.8 (23) 



Kelner et al. I (l2006l) have recently presented a new analytical 
expression (a modified exponential cutoff with exponent /I) for 
secondary gamma-ray spectra from inelastic proton-proton in- 
teractions based on Monte-Carlo simulations. Table |4]provides 
the results of fitting this function, in addition to several other 
functional forms to the data. A pure power-law model is clearly 
ruled out, the alternative spectral models provide significantly 
better descriptions of the data. For the modified exponential- 
cutoff shape, the exponent /5 was fixed at 1.0, 0.5 and 0.45. 
The latter yS value was chosen because it yields the smallest 
value. The cutoff energy of the gamma-ray spectrum was 
found to vary depending on values of ~ 18 TeV, ~ 3.7 TeV, 
and ~ 2.3 TeV were fit, and none of the three yS values is sta- 
tistically favoured over the other taking the value of as mea- 
sure. Under the assumption that indeed the VHE gamma rays 
are due to cosmic-ray protons interacting with ambient matter 
and subsequent tt" decay, one might get an idea of the param- 
et ers of the parent p roton spectrum. Following the approach of 
;lner et alj (120061) . proton cutoff values in the range of 50 to 
150 TeV with spectral indices a ranging from 1.7 to 2.1 would 
be fully compatible with the gamma-ray data presented here. 
Note, however, that a proton cutoff in the 100 TeV range does 
not mean the spectrum terminates at this energy. Especially in 
case of a hard power-law index a < 2 there would be a suffi- 
cient number of protons beyond the cutoff energy. 

Combining the data of three years it is possible to extend 
the gamma-ray spectrum up to energies beyond 30 TeV. Taking 
all events with energies above 30 TeV, the cumulative signif- 
icance is 4.8 cr. Table |5] lists all the flux points together with 
bin-by-bin event statistics. 



6. Summary & Discussion 

The complete H.E.S.S. data set of the SNR RX J1713.7-3946 
recorded from 2003 to 2005 is presented here. When analysing 
the data of different years separately and comparing them 
to each other, a very good agreement is found for both the 
gamma-ray morphology and the differential energy spectra. 
The H.E.S.S. telescope system obviously operates stably over 



the course of three years, if one takes known aging effects into 
account. 

A combined gamma-ray image using ~ 63 hours of 
H.E.S.S. observations in 2004 and 2005 was generated achiev- 
ing an unprecedented angular resolution of 0.06°. The mor- 
phology of RX J1713. 7-3 946 in VHE gam ma rays confirms 
its earlier characterisation ( Aharonian et al.| l2006b) of a thick, 
almost circular shell structure with the brightest regions in the 
northwest, very similar to the X-ray image of this source. The 
gamma-ray spectrum of the combined H.E.S.S. data over three 
years on RX J1713. 7-3946 extends over three orders of mag- 
nitude in energy. Although at the edge of sufficient statistical 
significance, the high-energy end of the gamma-ray spectrum 
approaches 100 TeV with significant emission (4.8cr) beyond 
30 TeV. Given the systematic uncertainties in the spectral de- 
termination at these highest energies and comparable statistical 
uncertainties despite the long exposure time, this measurement 
is presumably close to what can be studied with the current 
generation of imaging atmospheric Cherenkov telescopes. 

From the largest measured gamma-ray energies one can es- 
timate the corresponding energy of the primary particles. If 
VHE gamma rays are produced via ;7r" decay following inelas- 
tic proton-proton interactions, gamma-ray energies of 30 TeV 
imply that primary protons are accelerated to 30 TeV/0.15 = 
200 TeV in the sheU of RX J1713.7-3946. On the other hand, if 
the gamma rays are due to Inverse Compton scattering of VHE 
electrons, accelerated in the shell, off Cosmic-Microwave- 
Background photons (neglecting the presumably small contri- 
butions from starlight and infrared photons), the electron en- 
ergies at the current epoch can be estimated in the Thompson 
regime as Eg ^ 20 ^JE^ TeV » 110 TeV. At these large en- 
ergies Klein-Nishina effects start to be important and reduce 
the maximum energy slightly such that ~ 100 TeV is a realistic 
estimate. 

If one considers the functional representations found for 
the fit of the gamma-ray spectru m of RX J17 13.7-3946 (c.f. 
Table |4|, the basic findings of Aharonian et al. ( 2006bh are 
confirmed with improved statistics and increased energy cov- 
erage: a pure power-law spectral shape is clearly ruled out, 
alternative models like a broken power-law, a power with 
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Table 5. Flux points including relevant event statistics are listed for the spectrum of the combined H.E.S.S. data set, shown in 
Fig. ID For all 28 bins, the energy, the number of signal and background counts (ON and OFF), the normalisation factor a, the 
statistical sign ificance cr, the gamma-ray flux and the energy range of the bin are given. The significance is calculated following 
Li & Mai ( Il983i) . For the final bin, as it has only marginally positive significance, we list both the actual flux point and the 2cr 
upper limit (which is drawn in Fig.|4]l. Note that the energy and flux values given here are corrected for the variation of optical 
efliciency, as discussed in the main text. 



# 


E (TeV) 


ON 


OFF 


a 


cr 


Flux (cm ^s ') 


Range (TeV) 


1 


0.33 


5890 


5134 


1.00 


7.2 


(2.73 +0.38) X 10"'" 


0.30 - 0.37 


2 


0.40 


5583 


4797 


1.00 


7.7 


(1.48 + 0.19) X 10-'" 


0.37 - 0.44 


3 


0.49 


4878 


4010 


0.97 


10.5 


(1.13+0.11) X 10"'" 


0.44 - 0.54 


4 


0.59 


4202 


3409 


0.94 


11.6 


(7.22 + 0.63) X 10"" 


0.54 - 0.65 


5 


0.71 


3900 


2941 


0.94 


14.2 


(5.20 + 0.37) X 10"" 


0.65 - 0.79 


6 


0.86 


3682 


2833 


0.97 


11.9 


(2.56 + 0.22) X 10"" 


0.79 - 0.95 


7 


1.04 


3881 


2643 


0.98 


16.1 


(2.17 + 0.14) X 10"" 


0.95 - 1.15 


8 


1.26 


3982 


2758 


0.97 


16.0 


(1.40 + 0.09) X 10"" 


1.15 - 1.39 


9 


1.53 


4076 


2661 


0.98 


17.9 


(1.06 + 0.06) X 10"" 


1.39- 1.69 


10 


1.85 


3873 


2603 


0.97 


17.0 


(6.71+0.40) X 10"'^ 


1.69-2.04 


11 


2.24 


3452 


2251 


0.98 


16.8 


(4.50 + 0.27) X 10"'^ 


2.04 - 2.47 


12 


2.71 


3215 


2113 


0.98 


15.9 


(2.97 + 0.19) X 10"'^ 


2.47 - 2.99 


13 


3.28 


3075 


2081 


0.98 


14.6 


(1.95+0.13) X 10"'^ 


2.99 - 3.63 


14 


3.98 


2915 


2057 


0.98 


12.9 


(1.24 + 0.10) X 10"'^ 


3.63-4.39 


15 


4.81 


2537 


1721 


0.98 


13.1 


(8.91+0.68) X 10"" 


4.39-5.31 


16 


5.82 


2183 


1555 


0.98 


10.8 


(5.18 + 0.48) X 10"'^ 


5.31-6.43 


17 


7.05 


1961 


1525 


0.98 


7.9 


(2.77 + 0.35) X 10"'^' 


6.43 - 7.79 


18 


8.53 


1507 


1208 


0.98 


6.2 


(1.49 + 0.24) X 10"" 


7.79 - 9.43 


19 


10.33 


1211 


881 


0.98 


7.6 


(1.27 + 0.17) X 10 


9.43 - 11.41 


20 


12.51 


881 


664 


0.99 


5.8 


(6.69 + 1.15) X 10"'" 


11.41 - 13.81 


21 


15.14 


652 


551 


0.99 


3.2 


(2.58 + 0.82) X 10""* 


13.81 - 16.72 


22 


18.32 


473 


364 


0.99 


4.0 


(2.19 + 0.55) X 10"'" 


16.72 - 20.24 


23 


22.18 


327 


260 


0.99 


2.9 


(1.10 + 0.38) X 10"'" 


20.24 - 24.50 


24 


26.85 


220 


153 


0.99 


3.6 


(8.82 + 2.46) X 10"'^ 


24.50 - 29.66 


25 


32.50 


182 


110 


0.99 


4.3 


(7.70 + 1.79) X 10"'' 


29.66-35.91 


26 


47.19 


227 


180 


0.99 


2.5 


(1.15+0.47) X 10"" 


35.91 -63.71 


27 


81.26 


51 


37 


0.99 


1.5 


(2.36+ 1.55) X 10"'* 


63.71 - 113.02 


28 


169.79 


14 


11 


1.00 


0.6 

Upper Limit 


(3.77+|39j X 10"'^ 
1.6x 10"'« 


113.02-293.82 



energy-dependent exponent, and a power law with exponen- 
tial cutoff describe the data significantly better Assuming 
an exponential-cutoff shape, a "slow" cutoff with exponent 
j6 = 0.5, as suggested by detailed Monte-Carlo simula- 
tions dKelner et al.l 120061) . yields a perfect description of the 
data, however, different values of /3 cannot be distinguished, 
but would rather require better event statistics at the highest 
energies. 

Given the good agreement of the results presented here with 
the previously published ones, our restrictive conclusions re- 
garding the nature of the parent particles re main unchanged 



to those outlined in Aharonian et al. (2006b). Both scenarios 



with a leptonic, or hadronic primary particle distribution are 
able to accommodate an exponential-cutoff shape with an in- 
dex of /3 ~ 0.5. However, if the mean magnetic field in the 
SNR region is indeed strongly amplified by the shock to val- 
ues well beyond typical interstellar fields, the hadronic na- 
ture of the observed gamma-ray emission would be difficult 
to concea l and this latter emissio n scenario would be clearly 
favoured dBerezhko & Volkll2006l) . 



With the deep H.E.S.S. observations of RX J1713.7-3946 
we approach now energies, at which attenuation due to pair 
production on the Galactic inter stellar radiation fi eld begins to 
affect the gamma-ray spectrum ( Zhang et al.ll2006;) . At the cur- 
rently measured maximum energy this effect is negligible, par- 
ticularly since RX J1713. 7-3946 is neither in the direction of 
the Galactic Center (more than 12° angular separation in line- 
of-sight), no r is it at the distance wher e the interstellar radiation 
field peaks dMoskalenko et al.l 1200 6). RX J1713.7-3946 wifl 
therefore presumably not be the astronomical source, where we 
will obtain a clear observational confirmation of the attenuation 
of gamma rays due to the interstellar radiation field. However, 
RX J1713. 7-3946 remains an exceptional SNR in respect of 
its VHE gamma-ray observability, being at present the remnant 
with the widest possible coverage along the electromagnetic 
spectrum. The H.E.S.S. measurement of significant gamma-ray 
emission beyond 30 TeV without indication of a termination of 
the high-energy spectrum provides proof of particle accelera- 
tion in the shell of RX J1713.7-3946 beyond lO'"* eV, up to 
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energies which start to approach the region of the cosmic-ray 
knee. 
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